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Tests of gravity on large-scales in the universe can be made using both imaging and spectroscopic 
. - - surveys. The former allow for measurements of weak lensing, galaxy clustering and cross-correlations 

such as the ISW effect. The latter probe galaxy dynamics through redshift space distortions. We use 
. a set of basic observables, namely lensing power spectra, galaxy-lensing and galaxy-velocity cross- 

spectra in multiple redshift bins (including their covariances), to estimate the ability of upcoming 
surveys to test gravity theories. We use a two-parameter description of gravity that allows for the 
f~| , Poisson equation and the ratio of metric potentials to depart from general relativity. We find that 

^ ' the combination of imaging and spectroscopic observables is essential in making robust tests of 

gravity theories. The range of scales and redshifts best probed by upcoming surveys is discussed. 
, We also compare our parametrization to others used in the literature, in particular the 7 parameter 

modification of the growth factor. 

PACS numbers: 98.80.Es, 98.62.Sb 
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CLc I. INTRODUCTION 

O 

%^ ' General relativity (GR) plus the Standard Model of particle physics can only account for about 4% of the energy 
C/3 , density inferred from observations. By introducing dark matter and dark energy, which account for the remaining 96% 
, of the total energy budget of the universe, cosmologists have been able to account for a wide range of observations, 

from the overall expansion of the universe to various measures of large scale structure ^ . 

The dark matter/dark energy scenario assumes the validity of GR at galactic and cosmological scales and introduces 
exotic components of matter and energy to account for observations. Since GR has not been tested independently 
on these scales, a natural alternative is that GR itself needs to be modified on large scales. This possibility, that 
modifications of the law of gravity on galactic and cosmological scales can replace dark matter and/or dark energy, 
, has become an area of active research in recent years. Attempts have been made to modify GR with a focus on 
• ' galactic Q or cosmological scales 0,0, The DGP model in which gravity lives in a 5-dimensional space-time, 
I can produce a late time acceleration of the universe. Adding a correction term /(i?) to the Einstcin-Hilbcrt action 
• 3 also allows late time acceleration of the universe to be realized. 

' In this paper we will focus on modified gravity (MG) theories that are designed as an alternative to dark energy 
• • , (DE) to produce the present day acceleration of the universe. In these models, such as DGP and f{R) models, gravity 
at late cosmic times and on large-scales departs from the predictions of GR. By design, successful MG models are 
difficult to distinguish from viable DE models using observations of the expansion history of the universe. However, in 
■ general they predict a different growth of perturbations which can be tested using observations of large-scale structure 

ci ; (Lss) p, Hi i M M M MM MM\iAMMM,M- 

LSS in MG theories can be more complicated to predict, but is also richer because different observables like lensing 
and galaxy clustering probe independent perturbed variables. This differs from conventional DE scenarios where the 
linear growth factor of the density field fixes all observables on sufficiently large-scales. Theories of LSS in these 
modified gravity models are still in their infancy. Most studies have focused on probes of a single growth factor 
with one or a few observables. Recent predictions for discriminatory power of different observables could be found in 

m HUE [13, [11. 

We study tests of gravity that can be made with a combination of imaging and spectroscopic surveys Our emphasis 
will be on model- independent constraints of MG enabled by combining different observables. Carrying out robust 
tests of MG in practice is challenging as in the absence of a fundamental theory, the modifications to gravity are often 
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parametrized by free functions, to be fine tuned and fixed by observations. Recently the Parametrized Post-Friedman 
approach has been suggested as an attempt to describe a variety of gravity theories [2^ . 

In f|TT]we describe the ingredients of our modehng - parametrization of the MG models f ^II Ap . observables used 
for forecasting f ^ill B[) and covariances between observables ( mi Cp . In mil] forecasts for upcoming imaging and 
spectroscopic surveys are presented. We conclude in ijlVI The derivation of formulas for some of the covariances used 
in the main text are given in the Appendix |X] and results for alternative MG parametrization in Appendix [Bj 



II. METHOD 



A. Parametrization of modifications to gravity 

We are concerned with sub-horizon scales that satisfy the quasi-static, Newtonian approximation. In this regime, 
two effective functions characterize the departure of modified gravity theories with scalar perturbations from general 
relativity. We neglect additional fields that are expected to play a role on small, nonlinear scales to drive the theory to 
GR. Scalar perturbations in a homogeneous and isotropic universe can be described in the Newtonian (longitudinal) 
gauge by two potentials $(i,x) and ^'(i,x) as follows 29] 

ds^ = - [1 + 2*(t, x)] dt^ + a\t) [1 - 2$(i, x)] [dx' + rixfdif] , (1) 

where a{t) is the scale factor and r{x) the comoving angular-diameter distance. Throughout this paper we assume 
that the universe is spatially flat, so that r(x) = x- In GR, neglecting sources of the anisotropic stress in the energy- 
momentum tensor, the relation $ = ^E" holds (2§|. Moreover, the curvature potential $ is related to the mass density 
distribution /5(a,x) = p(a)(S(a,x) through the Poisson equation, which can be altered in MG theories. We assume 
that the Fourier-space analogue of the Poisson equation becomes 

-/c2$(a, k) = 'iT:a^Gg{k)p5{a, k), (2) 

where Gg{k) is the effective gravitational constant. The relation between the curvature potential <& and the Newtonian 
potential \1/ in MG theories is parametrized as 

m = |- (3) 



There are a variety of 
number of authors 
Friedman approach 
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parametrization of MG in the literature; the one described above has been suggested by a 



|37| . |39|. A more general description of modified gravity models via a Parametrized-Post- 
22l | extends to the superhorizon regime, but we opt for the two free parameters described above 
as we restrict ourselves to the quasi-static Newtonian regime. It is expected to be a generic feature of MG that the 
linear growth of structure becomes scale dependent and its time evolution gets changed with respect to the GR case 

The simplest possibility is that the functions g{k) and rj[k) do not depend either on time or on scale, hence we 
denote g{k) = go and ry(fc) = 770- In GR both parameters are unity, which is their fiducial value throughout the paper. 
Other parametrization will be discussed below; in some cases a specific scale or time dependence allows for easier 
constraints on modified gravity parameters. 

The growth of structure in a CDM dominated universe with MG is given by 

S(a, k) + 2H{a)5ia, k) - ^4TrGpS(a, k) = 0, (4) 
r]{k) 

where the expansion history, given by the Hubble parameter H = a/a, will be taken to be identical to that in the 
standard ACDM cosmology. The growth of structure, described by the evolution of ^(a, k) as given in Eqn. 21 
depends on the ratio of g{k) and ?7(fc). The linear growth factor D{a) = d{a,k)/6{ai,k), which describes evolution 
of matter density perturbations (see Eqn. relative to their initial values a = ai, is sensitive to this ratio. The 
relevant growth factor for matter peculiar velocities (in fact, its divergence) relative to the density evolution is given 

We also consider a popular parametrization of MG which is based on the growth exponent 7 in the growth rate 
function /(a) [IB, For GR with a cosmological constant it can be expressed as /(a) = flmia)'* , where 7 = 0.55 
and Jlm(a) is the total matter density parameter. For the DGP model the growth exponent is 7 = 0.68 [l6l|. Other 
parametrization are discussed in the Appendix. 
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FIG. 1: Sensitivity of the growth function D(a) and the growth rate /(a) to the modified gravity parameter (left panel), 
growth exponent 7 (middle panel), and the dark energy equation of state parameters uiq and w\ (right panel) as a function of 
redshift. 

The sensitivity of the growth factor D(a) and the growth rate function /(a) to 770 is shown in the Fig. [TJ The effect 
of MG appears to be significant as for the redshift range we are interested in (z < 1.5) the change of D{a) is ^ 2% 
and /(a) is ~ 0.5% if go is changed by 1% (the effect of 770 has the same effect but with opposite sign). If the growth 
exponent parametrization is employed the change in D{a) and /(a) is significantly smaller. For comparison we also 
show the response of D{a) and f{a) to the change of the dark energy equation of state parameters wq and wi defined 
a.s w — wq + wi{l a). 

B. Observables 

We consider the signal from weak gravitational lensing together with galaxy clustering in redshift space. In metric 
theories of gravity the light deflection angle a is given by the transverse gradient of the sum of the metric potentials: 
a{9) = V_l($ -I- (see e.g. [13] )• Therefore observed shapes of galaxies and their correlations as described by weak 
lensing power spectra are dependent on both metric potentials. The clustering of matter is governed by the growth 
equation ^ which is dependent only upon the Newtonian potential. Different combinations of the information from 
weak lensing and redshift space galaxy clustering has been used to forecast tests of MG models [l^, 113, . 

For weak lensing we use two observables - the correlations between shapes of galaxies (cosmic shear) quantified by its 
power spectrum Ckk,{1) and correlations between the foreground galaxy distribution and shapes of background galaxies 
(galaxy-galaxy lensing) described by the Cgi^lJ) cross-power spectrum. We use the convergence field k for simplicity 
as the power spectra defined below for k are identical to the shear power spectra [ll|. It is given by: k{9) = ^Vg a{6). 
The relation between the convergence and the metric potentials is given by the line-of-sight projection: 



<0) 



dz r{z)r{zs,z) ^ 
H(z) r(Zs) 



(5) 



where r(z) is the comoving angular-diameter distance between observer and a lens at redshift z. We take the sources 
to lie at redshift Zs- 

The metric potentials are related to the mass distribution as given by Eqns. ^ and ([3]), so the lensing power 
spectra can be expressed in terms of the three-dimensional mass power spectrum P^{a,k). In the small-sky-patch 
limit (we work with scales smaller than ^ 6° ot I > 30), the Limber approximation [54| gives 



where the lensing weight function 
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depends on the geometry and the redshift distribution of lensed galaxies dub/dz. The three-dimensional wavenumber 
k is given by A: = l/r{z). We use lensing tomography (531] by dividing the galaxy distribution into bins in redshift. 
Hence, instead of one projected power spectrum we obtain Nz{Nz + l)/2 power spectra CniKjil) that carry additional 
information about the growth of structure. Similarly for the galaxy-shear power spectra we have 
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dz, 
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Pssik,Zi)WL{z,Zj). 



(8) 



We assume that the distribution of galaxies Sg is a biased tracer of the mass distribution S but their relation is local 
and given by a bias factor b which may depend on time, 6g{zi) = b(zi)d{zi). A non-zero Cg-^j is obtained when 
galaxies gi are in front of the source galaxies, which requires i < j. We also compute the galaxy-galaxy projected 
power spectrum which will not be used as an observable but is required in making forecasts for MG by means of the 
Fisher matrix approach. It is given by 
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b\z,) 
r^{z,) 



H{zi) 



dzi 
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(9) 



where we have assumed that galaxies in two redshift bins are not correlated with each other, a good approximation 
for wide enough redshift bins. 

Note that the 'observed' power spectra differ from the spectra in Eqns. (O, ([5]), (H]) because the effect of discrete 
sampling of the underlying convergence and galaxy density fields should be taken into account. It leads to shot (shape) 
noise terms in the case of the 'observed' galaxy (shear) power spectra 



(10) 

(11) 

(12) 



where rig is the projected density of galaxies. The shape noise term, proportional to Ue, accounts for the intrinsic 

ellipticities of galaxies and its value (per component) is taken to be 0.4/\/2. The cross-power spectra Cg.,^- are immune 
to the discrete sampling noise. 

Modifications to GR enter the projected power spectra through the lensing-specific factor 5(fc)-^^^y^ which is 
responsible for the relation between the mass distribution and the metric potentials. Moreover, the evolution of 
structure as expressed by growth functions D{a) and /(a) is affected by changes in gravity as a result of modifications 
to Eqns. dH) and dS]). 

In physical space, the power spectrum of the galaxy distribution Pgg{k) is expected to be isotropic, but in redshift 
space peculiar velocities distort the distribution of galaxies along the line of sight. The radial component of peculiar 
velocities cause the observable redshift-space power spectrum Pgg\k, fik) to be 'squashed' along the line of sight on 
large scales (in the linear regime) and to produce pronounced 'finger-of-God' features on small scales (in the nonlinear 
regime) [H, [H^. The directional dependence of Pgg^ is given by fj,k = k^\/k, which depends on the angle between a 
wave vector k and the line-of-sight direction. 

Although the picture is more complicated in reality (see [s^ for a detailed discussion) , it is a good approximation to 
decompose the redshift space power spectrum in terms of three isotropic power spectra relating the galaxy overdensity 
Sg and peculiar velocities v: the galaxy power spectrum Pgg(k), the velocity power spectrum Pyy{k) and the cross 
power spectrum Pgv{k) as follows [52„ ,55| 



P^gHk, = [Pgg{k) + 2iilPg,{k) + 4P..{k)] Fik^ixlal 



(13) 



where the term F{k'^ii1al) describes non-linear velocity dispersion effects. We set F « 1, which is a valid on sufficiently 
large scales for forecasting purposes. As before, we assume that galaxies are biased tracers of mass and the bias is 
time-dependent but scale independent. What we refer to as the velocity is actually the velocity divergence, which is 
related to the mass distribution through the continuity equation in the linear regime S + V ■ v/a = (sec 58]). Even 
if gravity is modified the continuity equation stays unchanged as long as there are no components interacting with 
matter. 

The redshift space power spectrum in the form (jl3p shows a distinctive pattern in its angular dependence which 
is important in obtaining the component power spectra from Pgl^ ; it has been used to measure Pgg{k), Pgy{k) and 
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FIG. 2: Examples of the shear-shear and galaxy-shear power spectra for the DES (left panel) and a Stage-IV survey (right 
panel). /,From top: the galaxy-shear cross power spectrum CgK with foreground galaxies at z = 0.4 and background galaxies at 
2 = 1, the shear-shear auto power spectrum Ckk- The shape noise contribution to the auto power spectrum for z = 1 is shown 
as well (dashed) . Note that the shape noise and galaxy shot noise contribute to the variance of the power spectra. 



Pvvik) power spectra from the 2dF and SDSS galaxy surveys [53, [53]. Moreover, this decomposition is immune to 
the scale dependence of both the galaxy bias and growth functions as long as the angular structure is preserved [13] . 

The galaxy-velocity power spectrum Pg^ which we would like to use in constraining MG models is not a direct 
observable. We construct an estimator of Pgy band-power spectra as follows 



-Pgv (^i 



1 



W^9.(M)n'^^(fc,M), 



(14) 



where the summation is carried out over modes in a spherical shell of radius fc^, and the weight function Wg„(/i) is 
given by 



= ^P2(/i) 



135 



(15) 



in terms of Legendre polynomials Pi (/i) of order I (see Appendix A for motivation of this expression) . The summation 
in Eqn. p4p is carried out over modes contained in spherical shells in Fourier space which satisfy the following 
condition ki — Aki/2 < \k\ < hi + Afci/2. The volume Vk of this shell can be approximated by Vk = AnkfAki. The 
fundamental volume can be expressed in terms of the survey volume V^^ as Vp — {2Tr)^ /Vg. The number of modes in 
the volume Vk is then given by Nk — Vk/Vp. We assume that the power spectrum Pgp (fc, /^) does not vary significantly 
in the shell as a function of |A;|. We also assume that the survey is large enough that the fundamental volume is much 
smaller than Vk- Pgvih) is taken as an observable in the Fisher matrix analysis. Note that the redshift space power 
spectrum does not get any additional modifying factors in MG (except for the growth functions), as the continuity 

(s) 

equation on which the relation between Pgg and Pg-u is based remains the same. 

In the present work we focus on the galaxy- velocity cross power spectrum Pgv{k) and the information it can deliver 
about the growth of structure and galaxy bias. The bias factor appears linearly in our observables, both in projected 
power spectra as well as in Pgy{k), and is degenerate with the growth functions. However, when information in Pgy{k) 
is combined with tomographic measurements of the weak lensing signal it should allow for breaking this degeneracy 



C. Fisher matrix analysis 



In order to forecast the minimal attainable errors on cosmological parameters we implement the Fisher information 
approach [H, |3| including parameters that describe modifications to gravity. We treat as observables the following 
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band-power spectra, as described in ijllBI shear-shear power spectra Ck^kj (^m), galaxy-shear power spectra Cg-,^. (Im) 

and galaxy-velocity redshift space power spectra Pgvikm)- The data vector of the projected spectra can be written 

as = ^Ci^-f^j{l„i),Cg^K,j{l7n),Pgv{zi,kn)^, whcrc pairs of tomography bin indices : j > i with i,j = l,...,Nz 

denote independent power spectra. Each of them is comprised of Ni band-powers. There are also redshift space 
power spectra with Nk band-powers each. Therefore, the total number of observables in the Fisher matrix analysis is 
Nz{Nz + 1) X Ni + Nz X Nk- The projected spectra and the redshift space ones are independent so we may add their 
Fisher matrices or merge them into one data vector as above. 

The observable power spectra depend on a set of parameters pi whose uncertainties we aim to forecast. They are 
the following: (with fl\ adjusted to maintain spatial flatness), initial power spectrum slope rig, normalization 
of the power spectrum at the epoch of last scattering A^(fco = 0.002/Mpc), and a bias parameter in each redshift 
bin b{zi). In addition we use a two parameter [go, 770) or one parameter (7) description of departures from GR (see 
mi Al for definitions). We do not change either the dark matter or baryon physical density when other parameters are 
varied to avoid adding extra information from the change of the matter power spectrum shape. 

The Fisher matrix, which measures the curvature of the likelihood function in parameter space around its maximum, 
can be expressed for Gaussian distributed observables D^, as [1^ 

^..=i:i^c„v->0„D„)§, (16, 

The marginalized 68%-level error on a parameter pi is then given by cr^{pi) = [-^^^] j where F~^^ is the inverse of 
the Fisher matrix. In order to proceed with forecasting we need to compute the covariance matrices in Eqn. (fT6|) . We 
assume that the band-powers constructed from the convergence k and galaxy density Sg are Gaussian, which restricts 
the validity of our analysis to relatively large scales. An analysis of non-Gaussian effects on the lensing power spectra 
may be found in {45^ . 

The Gaussianity assumption allows us to express the covariance matrices of the 'observed' power spectra as follows 



Cov[C,^«^(/),a„,„(Oj = [C.,.Al)C.^.^{l) + C^,^^{l)C^^^,Sl)\ //skyA^(0, 

COV[C'3^«^(/),C'3„«„(01 = \Ca,g„Xl)Cn,nAOSrrr.+Cg,^^{l)Cg„^,,^{l)\ //skyiV(0 



Gov 



(17) 
(18) 
(19) 



where /sky is fraction of the sky covered by the survey and N{1) — ^;"'"'(2^ -I- 1) is the number of independent modes 
in a passband between ^min and ^max- We assume that modes in the power spectra are uncorrelated with each other. 

In our analysis the binning in multipoles I is logarithmic - we assume 15 bins in the range of multipoles from 
30 < Z < 1000. We choose this lower limit in order to be able to apply the Limber approximation when computing 
lensing power spectra. The upper limit is chosen to limit the nonlinear contributions to the power spectrum. Although 
at the high-Z considered there is a nonlinear enhancement to the power spectrum, we use the linear contribution as a 
conservative choice since both MG effects and biasing can become complex on small scales. We will show some results 
that include the nonlinear enhancement. 

In Fig. [2] we show example shear-shear and galaxy-shear power spectra along with the relevant errors from Eqns. 
P?)) and (HH]). The shape noise contribution to the observed spectra is also shown. The power spectra are flat as we 

— 1/2 

take only the linear evolution of mass perturbations into account. The errors on the power spectra scale with /g^y 
(see table HJ. 

We are interested in both Pgv{ki) and its covariance Cov[Pgy(ki), Pgv{kj)]. In order to compute the covariance 
one could use a standard approach and express the covariance by means of the survey's effective volume VcS as 
(ki) — 2V^^ Pg^{ki) [3^. Issl. [6]| . The effective volume accounts for the survey window function and in the sample 

variance limit approaches the physical survey volume Vs fs^ . [ssj . A drawback of this formula is that it underestimates 
the expected noise of the estimator Pgv{ki) and does not account for the effect of discrete sampling of redshift space 
with mass tracers (we thank R. Scoccimarro for drawing our attention to these issues). Let us consider the covariance 
matrix for the galaxy- velocity band-power spectrum Pgv{ki 



defined as 

COv[_Pgi, (fci) , i^gj, (fcj)] = {Pgi,[ki)Pgy[kj)) P gy (k^) P g^ (k 



Using the Gaussian assumption this is given by 

GOY[Pg^{ki),Pg,{kj)] = 



2(27r)3 



2TTkfAkiVs 



(20) 



(21) 
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FIG. 3: The galaxy-galaxy Pgg{k) (left panel), galaxy- velocity Pgvik) (central panel) and velocity- velocity Pvv{k) (right panel) 
power spectra at redshift 0.5 < z < 0.7 for the BOSS-I survey. Their expected statistical errors are shown in red. We assume 
15 bins per decade in k, and do not include nonlinear effects. 



Next we plug Eqn. (jl3p in to Eqn. (|2ip and integrate out the angular dependence to express the covariance matrix 
in terms of the component power spectra Pgg , Pgy , Pw as follows 



COv[Pgy{h),Pgy{kj)] = 



2(27r)3 



105 



2TTkfAkiV, ^ 18304 



3003 ( Pggihj 



+ 4420 Pggih) + — Pg4h) + 2940P„"„(fc,) 



+ 10P™(fc.) <; 147 ( Pggih) + ) + 226P3„(fc,) 



24185 
51 



(22) 
(23) 



The finite volume of the survey Vs is accounted for as well as the shot noise on small scales, which is inversely 
proportional to the mass tracer's density fig. A detailed analysis of the bias and covariance of band-power spectra 
estimators Pgg, Pgv and Pw can be found in the Appendix. Formulas for the covariance matrices of these power 
spectra are presented in Eqns. (|A16|) . (|A17|) and (jA18|) . 

In Fig. [3] we show galaxy and velocity power spectra, Pgg{k) and Pyy{k), and the galaxy- velocity cross power 
spectrum Pgy(k) together with their expected errors computed as described above. The plot shows as an example 
the power spectra from the BOSS-I survey redshift slice centered aX z = 0.6 and Az = 0.2. The errors on Pgy{k) 
are larger by a factor ~ 4 from these on Pgg{k). The good news about the cross power spectrum is that it does not 
suffer from the shot noise which dominates the galaxy power spectrum for k > 0.2/i/Mpc. On the other hand it is 
affected by sample variance; for k < 0.03/i/Mpc, Pgvik) has limited information. From Fig. [3] we see that the scale 
dependence of the fractional errors oi Pgvik) and Pwik) power spectra bears a typical u-shape with a minimum about 
k ~ 0.3ft-/Mpc. When making error forecasts we limit our calculations to 0.015 /i/Mpc < k < 0.15 h/Mpc to stay 
within the linear regime, and use 15 logarithmic band-powers in this wavcnumber range. 

We compare uncertainties on redshift-space power spectra to the results of [34( |. whose work is closely related to 
ours. At the fixed scale k = 0.05/i/Mpc our method yields errors on Pgg which are ~ 3 times smaller than presented 
in [33|. On the other hand errors on Pw in our analysis are ~ 2.5 larger than those in f33|. 

The total signal-to-noise ratio (S/N) for Pgvik) is presented in Fig. [4] as a function of redshift. We also show the 
S/N for the galaxy-shear power spectra Cg^ for the DES and a Stage-IV surveys assuming that source galaxies are 
at fixed redshift - Zs = 1.1 for DES and Zs = 1.9 for the Stage-IV. As expected, the S/N is the highest when lensing 
galaxies are about half way to the source galaxies. The most robust constraints on MG are expected when an imaging 
and spectroscopic survey both have high S/N over a common redshift range. For DES and BOSS-I this occurs at 
redshifts between about 0.4-0.5. 

It is worth noting that the S/N for Pgvik) is inversely related to the bias b of the galaxy sample. A joint analysis 
by combining different population galaxy samples may be helpful in beating down the sample variance on large scales 
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FIG. 4: Left panel: Redshift distribution of galaxies for imaging surveys specified in Table U] The distributions are normalized 
so that the maximum value is unity. Right panel: The total signal-to-noise ratio for the galaxy-shear power spectra Cqk. as a 
function of redshift of lensing galaxies. The background galaxies are located in the furthest bin - z = 1.2 for DES, z = 2.0 for 
a Stage-IV survey. The galaxy-velocity power spectrum signal-to-noise is also shown for a BOSS-I-type survey (extended to 
z = l.Q with the same sky coverage and galaxy density as the dotted curve). 
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TABLE I: Parameters of imaging surveys: sky coverage /sky in sq. degs., galaxy surface density rig per sq. arcmin., the zo 
parameter of the galaxy redshift distribution, and its mean {z). 



D. Surveys 

One of the main scientific goals of upcoming multi-color imaging surveys is to measure cosmological weak gravita- 
tional lensing. We consider two surveys of this kind: the Dark Energy Survey (DES) 48] which is expected to begin 
data acquisition in 2011, and a generic Stage-IV survey ^42,] whose example is the LSST survey [ioj . 

The surveys are characterized by sky coverage /skyj surface density of lensed galaxies and the galaxy redshift 
distribution. The sky coverage for the DES is taken to be 5000 sq. degs. and for the Stage-IV survey 20000 sq. degs. 

The redshift distribution of galaxies in the imaging surveys is assumed to have the form dn/dz oc z^exp — (z/zq)*^^^ , 

where the values of Zq for surveys under consideration are given in the table |T] and the distributions are shown in the 
Fig. H 

is) 

In order to measure the redshift-space power spectrum Pgg we consider spectroscopic surveys. The Baryon Os- 
cillation Spectroscopic Survey (BOSS) [51] will target Luminous Red Galaxies (LRG) up to redshift z ~ 0.7 and 
will cover a quarter of the sky. It will obtain spectra of 1.5 x 10^ galaxies, which results in the number density 
rig = 1.1 X 10~'*Mpc~'^. Based on the LRG sample from the SDSS survey one expects these objects to be biased by 
a factor of 6 ~ 2 with respect to the mass distribution. We assume that galaxies are uniformly distributed across 
the redshift range. In addition to the BOSS survey (called BOSS-I throughout the paper) we consider a futuristic 
version (dubbed BOSS-II here) with double the sky coverage compared to BOSS-I (to keep up with the sky coverage 
of Stage-IV survey), the same galaxy number density, and extending to redshift z — 1.1. 



III. RESULTS 



We quantify constraints on modifications to gravity by combined measurements from upcoming imaging and spec- 
troscopic surveys. Imaging surveys (see Table H]) provide us with the lensing signal via Ckk(0 and and CgK,{l). The 
lensing power spectra depend on the modified gravity parameters and {CgK{l) only) on redshift dependent galaxy 
bias, as given by Eqn. ([S]). The redshift space power spectrum Pgv(k) can be obtained from spectroscopic surveys. It 
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/sky 


•^min 


•^max 




ng 




BOSS-I 10000 


0.1 


0.7 


15.5 


1.1 X 10"" 


0.05 


BOSS-II 20000 


0.1 


1.1 


90 


1.1 X 10"'* 


0.14 



TABLE II: Parameters of the spectroscopic surveys: sky coverage /sky in sq. degs., minimum and maximum redshift limits 
of the survey, Zmin and Zmax, survey comoving volume in Gpc'^, mean spatial density of galaxies fig per Mpc'' and their mean 
projected density n^'^ per sq. arcmin. 




0.24 0.26 0.28 0.24 0.26 0.28 



FIG. 5: Forecast constraints on the effective gravitational constant go, the ratio of metric potentials 770, and flm for the 
DES imaging survey and the BOSS-I redshift survey. Forecasts based on the linear lensing power spectra (shear-shear and 
galaxy-shear) for DES are shown in light gray, and combined with the BOSS-I galaxy-velocity power spectrum in dark gray 
(blue). Light (red) inner contours show the forecasts if the nonlinear lensing power spectrum is used. Other parameters are 
marginalized over, and all contours show the 68% confidence level. 

depends on the ratio of go and 770 and on galaxy bias. 

Our fiducial cosmological model is given by = 0.26, the present value of the Hubble constant Hq = h x 
lOOkm/s/Mpc = 72km/s/Mpc, logarithmic slope of the initial matter power spectrum n — 0.96 and its amplitude 
A^(fco — 0.0G2/Mpc) = 2.41 x 10~^ [25|. The redshift dependent galaxy bias is a free parameter in each redshift 
bin with a fiducial value 6 = 2, corresponding to the Luminous Red Galaxy (LRG) sample from the SDSS. Fiducial 
values for parameters which describe modifications to gravity are set to their values in GR: go ^ rjo = I and 7 = 0.55. 
We use statistical priors on fim, ^mh^, ^bh^, and the power spectrum parameters and n as expected from the 
Cosmic Microwave Background measurements by the Planck satellite [1^. We also assume that the distance- redshift 
relation is unchanged from the standard ACDM cosmology. Therefore, MG enters the equations through the growth 
of structure and influences the weak lensing and redshift space power spectra, but does not affect the distance-redshift 
relation. Last but not least, we present uncertainties on cosmological parameters after uncertainties in the other 
parameters were marginalized out. 

In our analysis we use tomographic measurements which require binning of the lensing and redshift-space power 
spectra in redshift intervals. We assume bins with width Az = 0.2. Thus for the lensing spectra from the DES survey 
we use 6 redshift bins, while for for the Stage-IV survey we use 10 bins. For spectroscopic surveys we use 3 redshift 
bins for BOSS-I and 5 for BOSS-II. The projected density of spectroscopic galaxies is smaller by a factor of a few 
tens than for the imaging surveys, as given in tables [11 and [Til This affects the errors on Cgg large but does not effect 
CgK which is one of our observables. 

Our results on MG parameters constraints for the DES and BOSS-I surveys are shown in Fig. [51 If we consider weak 
lensing observables only (with CMB priors as described above) the one-sigma error on go and rjo is cr(.go) = 0.94 and 
cr(?7o) = 1.2. The weak constraints are due to the strong covariance between these parameters and redshift dependent 
bias. However, when combined with BOSS-I the constrains on go and r/o improve to cr((?o) = 0-12 and cr(?7o) = 0.15, as 
shown in the Fig. [51 Thus redshift-space clustering data enables us to beat down errors on MG parameters by about 
a factor of 8 (note that the constraints on flm come mainly from the CMB prior). The BOSS-I survey as presently 
planned will obtain redshifts for objects with z < 0.7. With the expanded redshift survey BOSS-II, the accuracy 
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FIG. 6: Forecasts as in Fig. [5] but for a Stage-IV imaging survey and the BOSS-II spectroscopic survey. 



would improve to a{go) — 0.086 and 17(770) = 0.11. 

In order to use information on small scales from lensing one needs to model the nonlinear evolution of the matter 
density. Nonlinear evolution boosts the lensing signal Ckk(0 by a factor of ~ 4 on scales I ^ 500 — 1000. Numerical 
simulations which provide nonlinear matter power spectra in MG theories arc in their infancy but the first attempts 
are encouraging, as shown by ^2d\ for the /(i?) models and [13, [4l| for the DGP model. However, there is no simple 
and model-independent parametrization of the nonlinear corrections to the matter power spectrum in MG models. 
This is due to the complexity of the evolution equations and the existence of additional fields that drive the theories 
to GR on small scales. 

Therefore we use only the linear power spectrum up to ^ = 1000. This underestimates the signal-to-noise in lensing. 
An alternative is to use the GR-based nonlinear power spectrum, as our fiducial model is GR and small deviations 
from the GR may not introduce substantial deviations in the nonlinear evolution. In Fig. [5] we also show predicted 
uncertainties in the MG parameters if we include such a nonlinear power spectrum in the modeling of the lensing 
spectra. The errors on both go and 770 drop by almost a factor of 2 to a-{go) = 0.069 and 17(770) = 0.078, respectively, 
compared to the linear case. This improvement is partly due to breaking the degeneracy between redshift dependent 
bias and MG parameters in the nonlinear regime. 

We also examine Stage- IV-type surveys like LSST, which will have greater depth and sky coverage. The predictions 
are presented in Fig. [51 Without spectroscopic information the constraints are (7(50) = 0.19 and 17(7/0) = 0.24. If 
we combine imaging data with the redshift-space power spectrum Pg„(fc) from the BOSS-II survey, the constraints 
improve to <y{go) = 0.048 and (7(770) — 0.062. Including nonlinear evolution in the lensing spectra leads to about a 
factor of two improvement. 

The MG parameters are strongly correlated with each other as shown in Fig. [T] This is mostly due to the 
dependence of the growth factor solely on the ratio of the MG parameters. The dependence of the lensing power 
spectra on a different combination of MG parameters (see Eqn. ([5])) is not a strong effect in practice; moreover, 
the change in growth can be compensated by a change in the redshift dependent galaxy bias which substantially 
enhances correlations. For a different MG parametrization like the one presented in the Appendix |B] one expects a 
weaker degeneracy. By combining lensing data with the redshift-space Pgvik) the MG parameters are more reliably 
determined, as their correlation with the galaxy bias parameters is significantly reduced (see Eqns. ^ and (HU). 

The importance of breaking the degeneracy between the galaxy bias and the growth of structure is highlighted 
in Fig. [7] by considering the dependence of uncertainties in go and 770 on the limiting projected scale Z^ax in the 
lensing power spectra and the limiting physical scale /cmax in the Pg^{k) power spectrum. For DES and BOSS-I 
we find that these constraints are insensitive to ^max if linear spectrum is used; the correlation coefficient between 
gt) and galaxy bias is r{go,bi) ^ 0.8. The degeneracy can be broken by using the nonlinear lensing power spectra 
which introduces scale dependent growth, improving the constraints by a factor of '--^ 2 and lowering the correlation 
coefficient to r(gQ, hi) 0.5. Another way to beat down errors is to increase the range of fc-modes in the redshift space 
power spectrum. However, extending the range of fc-modes to fcmax > 0.15 requires modeling nonlinear evolution and 
velocity dispersion effects. Throughout the paper we use linear lensing power spectra with Z^ax = 1000 and linear 
redshift space power spectrum with /cmax — 0.15 h/Mpc. In this case r{go, h) ~ 0.8. 
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FIG. 7: Correlation between modified gravity parameters for the DES and BOSS-I surveys (left panel), and for Stage-IV 
surveys (right panel). The error contours are as in Fig. [5] 
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FIG. 8: Achievable errors in the growth parameter 7 around a fiducial value 7 — 0.55. There are shown 68% confidence level 
contours for survey configurations as described in Figs. [5] and [6] There are shown 68% confidence level contours as expected 
from the DES and Stage-IV imaging surveys (gray) and when combined with spectroscopic surveys (blue). The contour (orange) 
shows predictions when the nonlinear matter power spectrum is accounted for. The braneworld DGP model has 7 = 0.68. 



The bias value of a given galaxy sample is related to the S/N ratio for the redshift space power spectrum, as 
discussed in mi CI Lower bias implies higher signal-to- noise for Pgy{k), which is counterbalanced by a lower S/N for 
CgK in the joint analysis. The net effect is that for the fiducial bias b = 1 the constraints on go and 770 are virtually 
unchanged compared to the b — 2 case which we have used in the analysis. If a less biased sample with b = 0.5 is 
chosen, the predicted errors increase by 3%. The fiducial bias value is therefore unimportant for the purpose of this 
work. 

We have used the {go,r]o) parametrization of MG so far. Another parametrization, based on the growth exponent 
(7), has been shown to be useful in distinguishing between gravity models 17|. Our results for 7 are presented in Fig. 
IHl It shows that using DES lensing data we expect to constrain 7 to (7(7) — 0.19, while Stage-IV survey lensing data 
achieve (7(7) = 0.07. When information about the redshift space power spectrum Pgy{k) is included the improvement 
is about a factor of two for DES combined with BOSS-I, and less than that for Stage IV surveys. Thus combining 
imaging and spectroscopic surveys is less useful for constraining 7 than the parametrization used in the rest of the 
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paper (compare Fig. 7 with Figs. 5 and 6). This is not surprising, as any measure of the growth of structure constrains 
this one-parameter modification. While it may indeed capture the relevant physics in some MG models, in general 
robust constraints require tests of both the Poisson equation and the ratio of metric potentials. These require both 
lensing and dynamical information, as illustrated in Figs. 5 and 6. 

The two models that have been extensively worked out in the literature are f{R) and DGP models; for these the 
ratio of potentials departs from unity by tens of percent in the quasi-static, Newtonian regime relevant to large-scale 
structure. The gravitational constant however is close to its value in GR if it is defined using the sum of metric 
potentials in the Poisson equation (rather than the usual definition, used in this paper as well, with the curvature 
potential) [22]. Thus rjo is the more sensitive parameter for testing MG if current models are used as a guide. 

IV. CONCLUSIONS 

We have used a Fisher matrix approach to test modified gravity models using imaging and spectroscopic galaxy 
surveys. The expansion history universe is not likely to be sufiicient to test MG models as it could be mimicked 
by an appropriately evolving dark energy equation of state. Here we use observable consequences of the evolution 
of perturbations to test gravity. In particular we use weak gravitational lensing measured from multi-color imaging 
surveys and dynamical information from spectroscopic surveys of galaxies. Lensing is sensitive to both metric po- 
tentials, whereas dynamical effects are driven by the Newtonian potential. Combining these probes provides robust 
tests of gravity, which is valuable given that we currently have a very limited number of specific models that are at 
all plausible. 

We use three simple parametrization of MG models and perform a joint analysis of shear-shear, galaxy-shear and 
galaxy-velocity power spectra. With a two-parameter description of modified gravity, we find that combining the 
three observables is essential to obtain strong constraints on gravity. We give predictions for Stage III and Stage IV 
surveys. We also compare the two-parameter description to the commonly used single parameter description (via the 
7 parameter), and to a scale dependent description used for braneworld models. 

Our results highlight the need for imaging and spectroscopic surveys to probe the same redshift range (for the 
lensing mass and galaxy distribution respectively) . Planned imaging surveys that reach redshifts of unity and beyond, 
and spectroscopic surveys that measure galaxy clustering at z ~ 0.2 — 0.6 are already well suited for probing modified 
gravity. With a careful selection of galaxy samples to compare cross-spectra, such surveys will allow us to perform 
robust tests of gravity. The question of how to use small scale information requires significant work, as different 
modified gravity models show a variety of nonlinear effects on scales below ^ 10 Mpc. It may be that models will 
need to be tested individually on these scales. Even so, measurements of the two parameters we have used would 
provide a consistency test of the GR plus a smooth dark energy scenario over a wide range of scales. 

We have focused on the most common MG parametrization with the gravitational constant and the ratio of the 
metric potentials being free parameters. These parameters turn out to be strongly correlated with each other and 
with the evolving, scale independent bias for the observables we have used. With external information on bias, or 
by including galaxy-galaxy spectra, this degeneracy could be broken (though the bias parametrization would need 
to be more complex as well with the inclusion of scale dependence). Alternatively, a different, physically- motivated 
parametrization may be better able to capture the dependence of observables on MG, such as the one discussed in 
the Appendix [B] ^33, [33] . The general properties of a useful MG parametrization, one that is able to get the most 
information out of a given set of observables, have been studied by :38i. 

Finally, in our analysis when we have included nonlinear evolution we have simply assumed it follows the predictions 
for GR. This provides one scenario for MG constraints; it may be optimistic as the degeneracy between MG parameters 
and scale independent galaxy bias gets lifted in this case. The existence of scale dependent bias on scales where 
nonlinear effects are important would lead to weaker constraints. On the other hand specific signatures of non- 
linearity would make it easier to distinguish models. Clearly more work is needed to include small scale information, 
realistic biasing schemes and additional observables such as galaxy power spectra, CMB lensing and the ISW effect. 

Acknowledgments 

We are grateful to Gary Bernstein, Alex Borisov, Mike Jarvis, Marcos Lima and Pcngjie Zhang for many useful 
discussions. We especially thank Roman Scoccimarro for discussions and help with the covariance calculations. This 
work is supported in part by NSF grant AST-0607667. MT is supported by World Premier International Research 
Center Initiative (WPI Initiative), MEXT, Japan, by Grand-in-Aid for Scientific Research on Priority Area No. 467 
"Probing Dark Energy through an Extremely Wide and Deep Survey with Subaru Telescope" and on young researchers 
(Nos. 17740129 and 20740119). 



13 



APPENDIX A: COVARIANCE MATRICES FOR THE REDSHIFT-SPACE POWER SPECTRA 

We consider distribution of mass tracers (galaxies) in the redshift space 6g^^ and its power spectrum Pgp which is 
defined as 

{6(;Hk)Sl^H'k')) = {2n)'P^lHk)SDik - k% (Al) 
where Soik) is the Dirac's delta function. In the discrete case (jAl[) becomes 

(5(^)(fc,)4^)(-fc,)> = i^p(^)(fc,)<5.„ (A2) 

where discrete Fourier modes Sg^\ki) have units of volume. We are interested in computing bias and covariance 
matrices for the galaxy-galaxy band-power spectrum Pgg{ki), galaxy-velocity Pgv{ki) and velocity-velocity Pyy{ki). 
First, let us define estimators of these power spectra 

i^xy(fc.) = ^5]M^xy(A^)Pi^)(fc,M), (A3) 

where X and Y stand for galaxy (g) or velocity (v) fields. The weight functions WxxifJ-) are given by 

5 27 

Wggi^i) = P0i^i) - -P2{P) + yP4(M), (A4) 

It: IOC 

Wg.{p^) = -P2{^l)~—P,{^^). (A5) 

31 S 

WM = — ^(m), (A6) 

where Pi (p) is the l-th order Legendre polynomial as a function of the azimuthal angle between a wavevector fc and 
a line of sight (/Xfc = fcy/A:). The expression (jA6|) agrees with the one found by |£0||. By ensemble averaging of (jA3|) 
one can show that weight functions (|A4p . (jASp and (|A6p together with (|A3p provide unbiased estimators for Pgg{ki), 
Pgv{ki) and Pwiki). In the continuous limit, if the fundamental volume Vp = {2Tr)^ /Vg is small compared to the total 
volume of the Fourier-space spherical shell Vk where averaging is carried on, we obtain 

{Pxvih)) ^^^^jdk 27rfc2 j dfl WxYif^) [Pggih] + 2fi^Pg„{k,) + /P™(fc,)] • (A7) 

Next, we use orthogonality relation for Legendre polynomials and after integration obtain 

{PxYih)) - -^^Pxy(fc,) = Pxvih). (A8) 

Therefore estimators (|A3P of power spectra Pgg , Pgv and Pw are unbiased. 

Now let us turn to computing covariance matrices for the introduced estimators. The formulation is general for all 
three band-power spectra we are interested in. The difference is in the weight functions. ^From the definition of the 
covariance matrix we have 

C0Y[PxY{h),PxYikj)] = {PxY{ki)PxY{kj)) - PxY{h)PxY{kj). (A9) 

We plug Eqn. (jASj) in Eqn. (|A9|) and obtain 

GoY[PxY{ki),PxY{k,)] + PxY{ki)PxY{kj) 

= ^Y.T.^^y^^')^^y^^'')^P9fik,^i)P^gf{k',^,'),) (Aio) 

where averaging is over all modes is spherical shells of radii k and fc' in the Fourier space. We can express the estimator 
of the redshift space power spectrum by means of (|A2p which leads to 

C0Y[PxY{ki),PxY{kj)] + PxY{h)PxY{kj) 

= 7^7^ E E Vl^xy(/^)W^xy(M')(4'^(fe)4^H-fe)4^)(fe')4^H-fc')). (AH) 
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Note that k denotes discrete Fourier modes contained in a spherical-shell region of the Fourier space for a given 
band-power fc^: ki — Aki/2 < \k\ < ki + Aki/2. We assume that the Sj^^k) is the Gaussian random field which allows 
to simplify (jAlip considerably by applying the Wick's theorem (see e.g. [53]). We obtain 



Cov[PxYih),PxYik,)] + PxY{h)PxY{k,) 

1 



J2 E WxY^WxYif/) (<5(^)(fc)<5(^)(-fe))(4^)(fe')4^H-fc')) 



N?. (27r)6 

^ ' k,fi k' ,fi' 

+ {sl''Hk)6^;Hk')){s(;\-k)6^;H-k')) + (4^) (fe)^^) (-&')) (4^)(-fe)4^Hfc')>] • (A12) 

By applying relation (|A2[) we obtain 

Cov[PxYih),PxYikj)] + PxY{h)PxY{k,) 

= ^ll^^y^^')Pafik,n) E WxY{p')P^,f{k',t^') 



+ ^2 E E WxY{^i)WxYi^^')[PtJik,^i)rs,^k'S,,-,> 



+ ]^EE^^^(^)^^^(^')[^i9^(^''^)]'^'«^^'W- (A13) 

^ k,fi k' 

In the second and third terms we introduced Kronecker delta-type symbol Sp^q which means that only pairs of modes 
which wavevectors are opposite contribute to the second term and only these which wavevectors are equal contribute 
to the third term. These delta functions make one summation in the second and third terms drop out. Moreover, all 
functions are even with respect to /z. Next, in the limit of continuous we derive 

Cov[PxY{h),PxY{k,)] - E^'^^'^^ W^ly(M)[^i^Hfc,M)]'. (A14) 



fe 



Finally, we obtain a general expression for the covariance of the Pgg{ki), Pgv{ki) or Pwih) provided validity of the 
decomposition (fT3|) 

Gov[PxY{h),PxY{k,)] = J^22l^y/ ^jl\ ^ w^yi^')P!;lHk,^^)■ (A15) 

By means of the relation (jAlSp it is straightforward to obtain desired expressions for the Pgg{ki), Pgvih) and Pwih) 
band-power-spectra covariance matrices which are presented below. Notice that we included the effect of discrete 
sampling of the mass tracer distribution by including the shot noise term l/fig. 
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COy[Pgg{h),Pg,{k,)] - 



COY[Pg,{k,),Pg,{k,)] = 



Cov[P„D(fci),-P™(fcj)] 
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APPENDIX B: ALTERNATIVE PARAMETRIZATION OF MODIFIED GRAVITY 

A useful modification to GR ■was proposed by l33], ■which seeks to describe changes in the potential-density rela- 
tionship on large scales "while leaving small scales unchanged from GR. The modification is expressed in the form of 
a power series in aH/k, "which is the ratio of the proper scale of perturbations a/k to the horizon size 1/H. The 
Fourier-space analogue of the Poisson equation is assumed to be modified as follo"ws 

-fc2$(a, k) = 4:7ra^Gg{k)p6{a, k), (Bl) 

■where g{k) = go{a) + gi{a)^. The relation (IBip converges to GR ■when go — I and for scales much smaller than 
the horizon size, i.e. aH/k ^ 1. Note that we consider only the first two elements of the power series. The linear 
term in the expansion (jBip is characteristic for brane-world inspired models like DGP. The linear term is absent in 
scalar-tensor models including f{R) models, where the first non-zero higher order term is quadratic in aH/k. For a 
thorough discussion of the parametrization and specific examples in different alternative gravity models see [3^ 1 . 
Similar to the curvature potential <&, the Newtonian potential can be modified as 

-kH{a,k) = 4:na^Gfi{k)pS{a,k), (B2) 

where = /io(a) + /ii(a)^. In GR g^ — fiQ ~ ^ a-nd the scale dependence of both potentials vanishes. Thus, 
we parametrize the departure from GR using 4 parameters go, gi, /xp and fii. They are used in the Fisher matrix 
analysis with fiducial values are ffo = Mo = 1 and gi = pi = 0. This type of modification is supported by the fact that 
the growth of structures is sourced by the Newtonian potential \E', so that po and contain information about the 
effect of MG on the observed matter distribution. The growth equation takes the form [l^l 

1,2 

(5(a, k) + 2H{a)5{a, k) + //(/c) — $(a, k) = 0. (B3) 

The results for MG parameters with Stage-IV and BOSS-II surveys are shown in Fig. [9l The imaging Stage-IV 
survey (with the usual CMB prior) could constrain the scale-independent part of the effective gravitational constants 
as (j{go) = 0.22 and cr(/io) — 0.075 (about a factor of 4 smaller than those achievable with DES.) If we add information 
from BOSS-II spectroscopic survey we obtain a{go) = 0.050 and cr(/io) = 0.018 (about 2.5 times smaller than for DES 
and BOSS-I. ) 
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FIG. 9: Achievable uncertainty on modified gravity parameters go, no and matter density ^Im for the future configuration of 
the Stage-IV imaging survey and the BOSS-II galaxy redshift survey. Predictions for the Stage-IV survey only are shown in 
gray, for the Stage-IV and BOSS-II combined - in blue. There are shown 68% confidence level contours. The dependence on 
the other parameters is marginalized out. 



In Fig. [TO] we present joint constraints on the scale independent and dependent terms in the modification. The 
uncertainties on the scale dependent terms are much larger (compared to the constant terms): (^{gi) = 3.52 and 
(t(/ii) = 1.39. The weak constraints on the scale dependent terms relates to the fact that this dependence is important 
on large scales approaching the horizon, where the signal-to-noise is small. The integrated Sachs- Wolfe effect would 
be more promising to test for modifications on the largest scales. 

On the other hand, even constraining the scale independent modification lets us distinguish a class of brane-world 
models like DGP gravity from the standard LCDM one. The values of parameters under examination for the DGP 
model are go = 1.25, gi = 0.5, fj,o — 0.75 /ii = 0.5 at z ~ 0.2 (and closer to GR values for higher redshifts ^32,]). 
Thus the DGP model is several a away from the GR (see [g^] for current constraints on DGP). The differences in the 
growth history in these two gravity models helps discriminate them (po is the most tightly constrained parameter). 

We have also explored time evolving parametrization of g and rj. With additional parameters, constraints generally 
get weaker, but if one chooses a specific fiducial time evolution (such as linear or quadratic in a) it can become easier 
to test MG models. We leave a detailed exploration of time and scale dependence for future work. See [3^ for a 
recent study using principal components in the scale and time dependence. 
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